Fermi gas with attractive potential and arbitrary spin in one-dimensional trap 
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A gas of ultracold ''Li atoms (eflective spin 1/2) confined to an elongated trap with one- 
dimensional properties is a candidate to display three different phases: (i) fermions bound in Cooper- 
pair-like states, (ii) unbound spin-polarized particles, and (iii) a mixed phase which is believed to 
have some resemblance to the FFLO pairing, ft is of great interest to extend these studies to 
fermionic atoms with higher spin, e.g., for neutral ''"K, *'^Ca, *^Sr or ^^'^Yb atoms. Within the 
grand-canonical ensemble we investigated the fi vs. H phase diagram for S = 3/2 {fi is the chemical 
potential and H the external magnetic field) for the ground state using the exact Bethe ansatz 
solution of the one-dimensional Fermi gas interacting with an attractive (5-function potential. There 
are four fundamental states: The particles can be either unpaired or clustered in bound states of 
two, three and four fermions. The rich phase diagram consists of these four states and various mixed 
phases in which combinations of the fundamental states coexist. Bound states of four fermions are 
not favorable in high magnetic fields, but always present if the field is low. Working within the 
grand-canonical ensemble has the following advantages: (1) A universal phase diagram is obtained 
by scaling with respect to the interaction strength and (2) possible scenarios for phase separation 
are explored within the local density approximation. The phase diagram for the superposition of a 
Zeeman and a quadrupolar splitting is also discussed. 

PACS numbers: 71.10.Pm, 36.40.Ei, 51.30.+i 



I. INTRODUCTION 

Spin-imbalanced ultracold ®Li atoms confined to dif- 
ferent geometries are spin-1/2 fermion systems display- 
ing the interplay of Cooper pairing and spin-polarization 
and have been the subject of several recent studiesji"— 
Confinement to nearly one-dimensional tubes can be 
achieved if the ultracold cloud of atoms is subjected to 
a two-dimensional optical lattice, which defines a two- 
dimensional array of tubesi^ The tubes can be regarded 
as isolated if the confinement by the laser beams is strong 
enough to suppress tunneling between tubes. The scat- 
tering between atoms under transverse harmonic con- 
finement is subject to a confinement-induced resonance.^ 
Fine-tuning this Feshbach-type resonance, the interac- 
tion between the fermions can be made attractive and 
its strength can be variedi^ The interaction is local and 
can be approximated by a (5-function potential in space. 
The confinement along the tube is roughly harmonic and 
weak; it can be locally incorporated into the chemical 
potential. Consequently, these systems of fermions are 
only locally homogeneous and within the local density 
approximation display phase separation with the varia- 
tion of the chemical potential along the tube.-^*? 

One-dimensional spin-1/2 gases have been extensively 
studied theoretically. M. Gaudin^ and C.N. Yang^° ex- 
tended Bethe's ansatz for the Heisenberg chain^^ and 
Licb and Liniger's results for the locally interacting gas 
of bosonsi^ to obtain the exact solution for a gas of spin- 
1/2 fermions interacting via a J-function potential. It 
was shown by Gaudin^ and later by Takahashii^ and 
Laii^ that for an attractive interaction in the ground 



state there are two classes solutions of the discrete Bethe 
ansatz equations, namely, real charge rapidities and 
paired complex conjugated rapidities. The former rep- 
resent spin polarized particles and the latter correspond 
to bound states of the Cooper type. There are then three 
possible homogeneous phases at very low T, the (1) fully 
spin-polarized state (only real charge rapidities), (2) a 
phase without polarization, where all particles are bound 
in Cooper-pairs, (only complex conjugated rapidities), 
and (3) a mixed phase in which unpaired spin-polarized 
particles coexist with Cooper pairs. The Cooper pairs 
are gapped (i.e. it requires a critical field to break-up 
the bound states) and display no long-range order. Simi- 
lar results were obtained for the Hubbard model with at- 
tractive C/ji^ii^ The mixed phase has been interpretedi^ 
as the one-dimensional analogue of the Fulde-Ferrell- 
Larkin-Ovchinnikov (FFLO) state.— For an isolated tube 
there is no long-range order of the pairs and hence no 
order parameter; however, the Cooper-pair correlation 
function acquires a phase in the mixed state, that is be- 
lieved to be reminiscent of the space modulation of the 
order parameter in higher dimensions. A coupling be- 
tween tubes increases the effective dimension of the sys- 
tem so that long-range order can arise and it is believed 
that this could lead to the realization of the FFLO phase 
in an ultracold gas of atoms^^ii? FFLO related phases 
have been observed in the strongly anisotropic heavy elec- 
tron compoundiS'i'^o CeCoIns (the interpretation is still 
controversial^^) and in the quasi-two-dimensional organic 
compounds A-(BETS)2FeCl4 and (TMTSF) 2 0104^^1^ 

Tubes with ultracold gases of atoms provide the unique 
possibility to study fermion systems with a spin larger 



than 1/2, e.g., '"'K (spin 9/2), "SCa (spin 7/2), «^Sr (spin 
9/2) or ^^■^Yb (spin 5/2) atoms. With an attractive in- 
teraction atoms with spin S can form bound states of 
up to (25 + 1) particles, extending this way the concept 
of Cooper pairs to larger clusters.^'* Consequently, the 
phase diagram will have more possible pure and mixed 
phases. In this paper we investigate the phases that can 
arise in the ground state using the Bethe ansatz solution 
of the one-dimensional fermion gas with (5-function po- 
tential. Sutherland^^ generalized M. Gaudin's^ and C.N. 
Yang's^° Bethe ansatz solution (for spin 1/2) to an arbi- 
trary number of colors iV = 25* + 1 [SU(iV)-symmetry]. 
The classification of states and the thermodynamics of 
the gas have been derived long ago by one of the present 
authors^fi and the ground state equations and elementary 
excitations have been previously studied in Refs. .27..,28j. 
for both attractive and repulsive potential and arbitrary 
number of colors. Introducing different chemical poten- 
tials for each of the colors, these results are valid for an 
arbitrary level splitting of the iV-fold multipletj ^^'^^i^^ 
The results for the i5-function potential model, as well as 
other integrable models, have been extensively reviewed 
by one of the present authorsi^ In this paper we use the 
solution to study the phase diagram for a gas of fermionic 
atoms constrained to a tube. To address a concrete exam- 
ple here we consider S =^ 3/2; the numerical effort and the 
complexity of the phase diagram rapidly increase with S. 
For S = 3/2 there are fours basic states, namely, bound 
states of four, three and two particles, and unbound par- 
ticles, and the corresponding mixed phases, which can 
have up to four coexisting basic states. As a consequence 
of Pauli's exclusion principle, the bound states must in- 
volve particles all with different spin-components; oth- 
erwise the wave function has a node and the J-function 
potential cannot be active. 

Recently, in the context of ultracold atoms, the Aus- 
tralian National University group^Sr^^ reconsidered the 
SU(A^) (5-function potential model with attractive inter- 
action. In Ref. [3^1, using the thermodynamic Bethe 
ansatz, the dressed energy potential are re-derived for 
"spin-1 fermions", and used to discuss the phase dia- 
gram for a superposition of a Zeeman and a quadrupo- 
lar splitting. In Ref. [SJl, these results are extended 
to spin-3/2 fermions for fixed density, strong attractive 
coupling (Tonks-Girardeau gas limit) and Zeeman as well 
as quadrupolar splittings. This case makes contact with 
the present paper. The ground state energy, critical fields 
and the low-T specific heat for the SU(iV) case are ob- 
tained in Ref. [S^j- The thermodynamics for the SU{?>) 
case has been studied in Ref. [33 . Most results of these 
four papers^"— are already derived in Refs. [26l - l29| . 

Note that the two-body interaction for spin larger than 
1/2 does not necessarily have to have SU(A^) symme- 
try. Integrable one-dimensional continuum models for 
the low-density limit displaying pairing have been con- 
structed for bosonic and fermionic systems. In Ref. [3J| a 
model for spin-1 bosons with exchange interaction is pro- 
posed and solved exactly via nested Bethe ansdtze for the 



ground state and thermodynamics. An extension of this 
model to SO{5) symmetry for spin-3/2 fermions has been 
proposed and solved in Ref. [35;]; the authors obtain the 
thermodynamic equations and discuss the spectrum of el- 
ementary excitations. Further extensions to models with 
hidden Sp{2s + 1) and SO {2s + 1) symmetries for high 
spin-s fermions and bosons, respectively, can be found in 
Ref. [36j. The influence of a pure quadratic Zeeman ef- 
fect (quadrupolar splitting) on the Mott-insulator phases 
of hard-core one-dimensional spin-3/2 fermions has been 
studied via DMRG, leading to a rich phase diagram. "^^ 

There are several other theoretical studies of ultracold 
spin-1/2 atoms in one-dimension. The direct imaging of 
the density profiles of the spatially modulated superfluid 
phases in atomic fermion systems were obtained by solv- 
ing the Bogoliubov-de Gennes equationj^ The pairing 
states were investigated on a lattice by means of the den- 
sity matrix renormalization group method in Ref. 3y]; 
this study leads to a fourth possible phase (in addition 
to the paired, unpaired polarized and the mixed phases) 
consisting of a metallic shell with free spin-down (i.e. 
reversed spins) fermions moving in a fully filled back- 
ground of spin-up fermions. The crossover from three- 
dimensional (FFLO phase) to one-dimensional (mixed 
phase) behavior is addressed in 140], where the phase 
diagram for a weakly interacting array of tubes is calcu- 
lated. A quantum Monte Carlo study of one-dimensional 
trapped fermions with attractive contact interactions was 
presented in [4l|. Finally, using the Bethe ansatz the 
low temperature thermodynamics was calculated in Refs. 
di] and m. 

The rest of the paper is organized as follows. In Sect. II 
we present the model and the discrete Bethe ansatz equa- 
tions for perodic and open boundary conditions. The 
elongated traps correspond to open boundaries rather 
than to the usually employed periodic boundary condi- 
tions. This difference gives rise to surface or boundary 
terms, which, however, are not relevant in the dressed 
energy potentials in the thermodynamic limit. In Sect. 
Ill we present the numerical solution of the Bethe ansatz 
equations, the phase diagram for a Zeeman splitting and 
the local density profile along the trap. In Sect. IV we in- 
vestigate the case where in addition to the Zeeman effect 
there is a quadrupolar splitting. Although it is not clear 
if nonlinear Zeeman splittings are of relevance to ulta- 
cold atoms in one dimension, it is an instructive situation 
to study which has been considered in Refs. |30ll34ll37 1 . 
Conclusions are presented in Sect. V. 



II. MODEL AND BETHE ANSATZ 

The Hamiltonian for a gas of nonrelativistic particles 
with (25'-|-l) colors (spin S) interacting via an attractive 
(5-function potential is 
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where Xi are the coordinates, Np is the total number of 
particles and c is the interaction strength. By fine-tuning 
the confinement-induced resonance^ the interaction can 
become attractive and its strength can be varied. Here 
h /2m, where m is the mass of the particles, has been 
equated to 1, or alternatively it has been scaled into Ti 
and c. 



Bethe equations for periodic boundary 
conditions 
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A^ = kj and L is the length of the box. The energy and 
the momentum of the state are given by 

E -Jlk^ . P = Y.^j ■ (6) 
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i=i 



The states of the coordinate Bethe ansatz are plane 
waves constructed from the two-particle scattering ma- 
trix. This scattering matrix satisfies the so-called Yang- 
Baxter triangular relation, which is a necessary condition 
for integrability. As a consequence of the triangular rela- 
tion many-particle scattering processes can be factorized 
into two-particle processes and the order in which the 
individual scattering processes take place can be inter- 
changed (the order becomes arbitrary). 

The generalization of the Gaudin-Yang^iiiS solution to 
more than two colors-- consists of an iterative application 
of the Bethe- Yang hypothesis (generalized Bethe ansatz), 
such that one color is eliminated at each step, leading to 
N = 25* -|- 1 nested Bethe ansdtze. Each Bethe ansatz 
gives rise to a new set of rapidities, {kj}, j = 1, ■ ■ ■ , Np 

for the charges (coordinate Bethe ansatz) and {Aq }, I = 
1, • • • , A^ — 1, with a = 1, • • • , Af '^'^ for the internal degrees 
of freedom (spin). Here Af '^'^ is the number of rapidities 
in the l*^ set and a is the running index within each set. 
All rapidities within a given set have to be different to 
ensure linearly independent solutions. Consider fermions 
of spin S with Zeeman splitting and let us denote by 
Ns-m the number of particles with spin component m. 
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such that Np > M^^) > • • ■ > M^'^^l > 0. As a conse- 
quence of the SU{N) invariance of the model the nested 
Bethe ansdtze for periodic boundary conditions yield 
the following sets of coupled equation o^^i^^'^^ 

exp{zk,L) = n e{k, - A^) , j = 1, ...,iVp (3) 
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B. Betlie equations for open boundary conditions 

Eqs. Q are derived for the standard periodic bound- 
ary conditions. Tubes, however, are not periodic and 
better represented by open or refiecting boundary con- 
ditions. A particle reaching the boundary is then re- 
fiected undergoing kj -^ —kj but without changing its 
energy. The corresponding reflection matrix satisfies re- 
flection equations with the two-particle scattering matrix, 
extending the Yang-Baxter equations. All matrices can 
be diagonalized simultaneously."'^''*^ The total length of a 
period is now 2L, where L is the length of the trap. The 
open boundary Bethe equations for the present model 
were derived previously by Oelkers et al^ It is conve- 
nient to write the Bethe equations in a form similar to 
Eqs. 11]) by letting the indices j and a run from —Np 
to Np and -M^') to M('), respectively41i^ The Bethe 
Ansatz equations for open boundary conditions are then 
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Hence, there are twice as many rapidities and the box is 
also twice as large, leaving the density of rapidities un- 
changed. The main difference between open and periodic 
boundary conditions are then the independent factors in 
Eqs. ([7]) and ©j which contribute with 1/L terms to 
the rapidity densities. This is very similar to the effect of 
magnetic impurities in a chain. Also, for periodic bound- 
ary conditions the Bethe states are plane waves, while for 
open boundary conditions they are standing waves. The 
energy and momentum are still given by Eq. ([5]). 



C. Classification of states and energy potentials 

For an attractive interaction and large L, the solutions 
of the discrete Bethe equations can be classified accord- 
ing to (i) real charge rapidities, belonging to the set {fc^}, 
associated with unpaired propagating spin-polarized par- 
ticles, (ii) complex spin and charge rapidities, which cor- 
respond to bound states of particles with different spin 
components, and (iii) strings of complex spin rapidities, 
which represent bound spin states. ^^'^^ States in class 
(iii) are not represented in the ground state; these states 
correspond to excited states and are not considered here. 
This classification of states is completely analogous to 
that of the Anderson impurity of arbitrary spin in the 
U -^ oo limit^^'^'^ and the one-dimensional degenerate 
supersymmetric t — J model/Si 

Since only particles with different spin components are 
scattered, i.e. experience an effective attractive interac- 
tion, we may build bound states of up to (25" -I- 1) parti- 
cles. A bound state of n (71 < A^ = 25* + 1) is character- 
ized by one real ^("~^' rapidity and in general complex 
A*-'-' rapidities, Z < n — 1, given by 



j^(i) ^^("-i)+ip|c|/2 , l<n 
p — — (n — / — 1), —{n — ? — 3), • 
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These spin and charge strings form classes (i) and (ii), 
which are present in the ground statei^ The real rapidi- 
ties ^("~i) have all to be different and satisfy the Fermi- 
Dirac statistics, i.e. the states are either occupied or 
empty. (For S* = 1/2 the bound states are frequently 
called Cooper pairs, although this analogy is not rig- 
orous.) In the ground state the rapidities are densely 
distributed in the interval [-~Bi,Bi] and we denote with 
^^'■' (C): ' — 0: I7 ■ ' ■ J 25 the dressed energy potentials (en- 
tering the Fermi-Dirac distribution). For S = 3/2 the 
four energy potentials satisfy the following coupled linear 
integral equations (for larger spin see Refs. [26,28.29] ): 
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and /i; is the chemical potential for the bound states in- 
volving [l + 1) particles. The jjli determine the integra- 
tion limits Bi through the condition that e'^^\-iiBi) = 0, 
since occupied states correspond to e'-'-'(C) < and for 
empty states the quantity is positive. For a pure Zee- 
man splitting the chemical potentials are /ig = fi + 3H/2, 
fii=fj, + H,ii2=tJ- + H/2 and fi^ — /i, where H is the 
Zeeman energy and /x the chemical potential for the total 
number of particles f2§ 

Note that if all the e*^'^ are rescaled to £'''/c^j M to 
n/c^, H to H/c^, all Bi to Bi/\c\ and £_ to ^/|c|, the 
equations are universal, i.e., independent of the magni- 
tude of \c\. Hence, within the framework of the grand 
canonical ensemble, without loss of generality, it is suffi- 
cient to present the results for \c\ = 1. The problem has 
then only two parameters, namely, H and fi. 

An ultra-cold atom system is inherently inhomoge- 
neous since the diameter of the tube gradually changes 
with position from the center of the trap to its bound- 
aries. As a consequence of the changing diameter of the 
tube, the quantization in the plane transversal to the 
tube gradually changes the zero of energy. This change 
can be represented by a harmonic potential, so that the 
actual local chemical potential /i is a function of x given 

by 



fi{x) + ^^mujf^gX^ = const 



(15) 



Within the local density approximation, it is ^{x) that 
enters the Bethe equations, pT|) - (ITil) . The solution is 



then exact for the one-dimensional system, but approxi- 
m.ative for the trap. This approximation^i^i^ is expected 
to be good since the variation of /x with x is sfow. The 
approximation neglects the quantization of the harmonic 
confinement, which is treated classically and locally in- 
corporated into the chemical potential. 

Fluctuations in the particle density arise due to the 
x-dependence of /i and due to possible weak Josephson 
tunneling between tubes. It is then necessary to solve the 
Bethe equations in the grand canonical ensemble rather 
than for the canonical ensemble, i.e. at constant number 
of particles. The results in the canonical ensemble may 
then no longer be universal, i.e. the scaling with |c| dis- 
cussed above is no longer valid, because of the constraint 
of fixed number of particles, which is the case treated in 
Ref. [31]. 



III. PHASE DIAGRAM FOR ZEEMAN 
SPLITTING 



In this section we discuss the phase diagram for a pure 
Zeeman splitting of the levels. The set of equations pTjl - 
(114)) is solved numerically by iteration. In this case, the 
energy potential £^^' corresponds to unpaired particles 
with spin-component Sz = 3/2; the energy e^^^^ to bound 
pairs with spin-components Sz — 3/2 and Sz = 1/2; the 
potential e^^-' to bound states of three particles of spin 
components 5*2 ~ 3/2, Sz = 1/2 and Sz = ^1/2, re- 
spectively; and finally e'^^ to bound states of four parti- 
cles all with different spin-components. We denote these 
states with Roman numbers, I, II, III, and IV, respec- 
tively. These states can coexist in mixed phases, for ex- 
ample we denote with I-I-IV the coexistence of unpaired 
and bound states of four particles and with I+II+III a 
phase where all states except four-particle bound states 
are present. 

The phase diagram of fi{x) vs. H for \c\ = 1 is shown 
in Fig. 1. Other values of |c| can be reduced to this phase 
diagram by adequately scaling /i and H. Note that all 
phase boundaries are given by the zero of some energy 
potential. The phase boundaries are then crossover lines, 
which are accompanied by a square-root singularity of 
one of the densities of states (one-dimensional van Hove 
singularity);^ For small magnetic fields all particles are 
bound in four-particle bound states (generalized Cooper 
pairs). The shaded area is the region where all bands 
are empty (system without particles). With increasing 
field other phases become realized. At very large mag- 
netic fields and/or for low values of fi (small number of 
particles) the phase IV is not favorable. For large ^ and 
intermediate magnetic fields all for bands are populated 
and hence spin-polarized unbound particles coexist with 
all possible bound states. This phase diagram can be 
compared to the one obtained in Ref. l3l| (Figs. 3 and 
4) for fixed number of particles (canonical ensemble). For 
a pure Zeeman splitting and as a function of field these 
authors obtain crossovers from phase IV to I-I-IV to I. 
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FIG. 1: Ground state phase diagram /i vs. H for a homoge- 
neous fermion gas of spin S = 3/2 with \c\ = 1. The shaded 
area corresponds to the empty system (no particles). The 
Roman numbers denote the number of particles involved in 
a bound states. Regions with more than one Roman number 
have coexisting or mixed phases with coexisting order. Note 
that in the vertical axis /x is a function of x as given by Eq. 
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FIG. 2: Expanded view of the ground state phase diagram 
/i vs. H of Fig. 1 for a homogeneous fermion gas of spin 
S = 3/2 with |c| = 1 to show the multiple phases for the low- 
density region. The shaded area corresponds to the empty 
system (no particles). 



A constant number of particles corresponds to a curve 
/i(-ff) in Fig. 1. The strong interaction limit considered 
in Ref. [31| corresponds to a low particle density. The 
sequence of phases we then expect from Fig. 1 is also IV 
to I-I-IV to I, in agreement with Ref. [3l|. 

Fig. 2 shows the expanded view of the region for small 
/i and 7J, which displays multiple crossovers. For larger 
/i there is a small region where the phase I-I-III-I-IV is 
stable. As mentioned above, the harmonic confinement 
of the trap can be treated quasi-classically and can be 
absorbed into the chemical potential via Eq. (|15p .^^ The 



chemical potential then decreases as we move from the 
center of the trap towards the boundaries. Hence, we 
move downward along a vertical line on the phase dia- 
gram. This can give rise to phase separation along the 
length of the trap. For instance, for H — 2 for a suffi- 
ciently high density of atoms, at the center of the trap 
the phase with all bound states coexisting (I-t-II+III+IV) 
would be favored, then moving towards the end-points (in 
either direction) of the trap first the four-particle bound 
states disappear (phase I-I-II-I-III), then the bound states 
of three particles (phase III) are depopulated and polar- 
ized unbound particles coexist with bound pairs (phase 
I+II), and finally a fully polarized gas phase (I) is possi- 
ble. For H = I, on. the other hand, we again could have 
the I+II+III+IV phase at the center of the trap and by 
moving to the boundaries we would observe the I-fll-l-IV 
mixed phase, then the I+IV phase and finally unbound 
fully polarized atoms (I). 
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FIG. 3: Density profile within the local density approximation 
for H = 1.5, /i(0) = -0.7 and /i(L/2) = -2.0. The position 
along the trap is given by Eq. p6[) . The three crossovers 
between phases are shown by the thin vertical lines. The 
densities n, of bound states of g + 1 particles (or polarized 
unbound particles if g = 0) are given by the solid (ns), dashed 
(n2), dotted (ni) and dash-dotted (no) curves. 

The local density profile as a function of x for the dif- 
ferent phases for H = 1.5 is displayed in Fig. 3, where 
Eq. ((T5|) was used to parametrize the chemical poten- 
tial in a trap of length L. Given /x(0) and ^{L/2)^ i.e. 
the chemical potential at the center and boundary of the 
trap, the position along the trap is given by (from Eq. 



m) 



xi{L/2) = ^/w^~lmJw^m^lm\ • m 

The density function of the rapidities is obtained from 
the dressed energies e''^ (^) by differentiation with respect 
to /i, i.e.^>S9 



p'^\0 



p'^KO 



1 de^^\0 
2tt dfi 



where p^''\(,) is the particle density and pjf (C) the cor- 
responding hole density for bound states involving q + 1 
particles. The integral equations satisfied by the density 
functions is similar to the one for the dressed energy po- 
tentials. After solving these equations numerically, the 
number of bound states (or polarized unbound particles 
if g = 0) per unit length is obtained from 



d^P^'HO 



(18) 



Note that the densities vanish with a square-root singu- 
larity that is characteristic of one-dimensional van Hove 
singularities as can be seen in Fig. 3. 



IV. PHASE DIAGRAM FOR ZEEMAN AND 
QUADRUPOLAR SPLITTING 



For S > 1/2 level splittings other than the Zeeman 
effect are possible. Here we consider a quadrupolar split- 
ting superimposed with the Zeeman effect. Nonlinear 
Zeeman splittings have been considered previorisly^2i^i22 
in similar contexts. In particular, in Refs. |3C)ll3l| Bethe 
ansatz equations for the same model but a fixed number 
of particles (canonical ensemble) were employed, while in 
the present paper we use the grand canonical ensemble 
(variable number of particles). For S = 3/2 the Zeeman 
and quadrupolar splittings are not the most general case, 
since also octupolar splittings are possible. It is not clear 
if such a situation (quadrupolar and octupolar splittings) 
can be realized experimentally for ultracold atoms in ID, 
but the problem is theoretically sufficiently interesting to 
be addressed. The Hamiltonian for the level splitting is 
given by 



Hspi = -HS, + D[3Sl - S{S + 1)] 



(19) 



where D can be either positive or negative. To be spe- 
cific we consider again the case S = 3/2. Note that the 
Hamiltonian ^ commutes with Tispi, so that the Bethe 
states also diagonalize H + Tispi ■ Keeping D fixed and as 
a function of H, 'Hapi displays two level crossings; hence, 
we need to consider three regions, namely. Region (i) 
{H < 311)1), Region (n) (3|D| <H< 6|D|) and Region 
(in) {Q\D\ < H). For Z? > the bound states I, II, HI 
and IV are then composed by particles with the following 
S:, components: In Region (i), (1/2), (1/2,-1/2), (1/2,- 
1/2,3/2), and (1/2,-1/2,3/2,-3/2), respectively; in Region 
(ii), (1/2), (1/2,3/2), (1/2,3/2,-1/2), and (1/2,3/2,-1/2,- 
3/2), respectively; and in Region (iii), (3/2), (3/2,1/2), 
(3/2,1/2,-1/2), and (3/2,1/2,-1/2,-3/2), respectively. Re- 
gion (iii) is then similar to the case of a pure Zeeman 
splitting. Hence, at the crossovers the character of the 
bound states changes. The corresponding chemical po- 
tentials ^i in Eqs. ([TT|) -([T4t are then 



(17) 



Region (i) {H < 3D) 

Ho = P + 'iD + 3i//2, Hi = fi 

H2= P + D + H/6, H3 ^ P 
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FIG. 4: Ground state phase diagram /i vs. H for a homo- 
geneous fermion gas of spin 5* = 3/2 with |c| = 1 and a 
quadrupolar splitting D\iSl - S{S + 1)] for D = 1/3. The 
shaded area corresponds to the empty system (no particles). 
Level crossings of the Zeeman and quadrupolar terms occur 
at if = 3D and H = 6D. 




FIG. 5: Ground state phase diagram fi vs. H for a homo- 
geneous fermion gas of spin S — 3/2 with |c| = 1 and a 
quadrupolar splitting D[35'f - S{S + 1)] for D = -1/3. The 
shaded area corresponds to the empty system (no particles). 
Level crossings of the Zeeman and quadrupolar terms occur 
at H ^ 3\D\ and H = 6|D|. 



Region (ii) {3D < H < 6D) 

^10 = ^i + 3D + 3H/2, ^1 = /i + iJ 

^2 = Ai + -D + H/6, A^s = A^ 



Region (iii) {6D < H) 

f^o = f^ + SD + 3H/2, fii = fi 

Ai2 = /i - £> + H/2, M3 = M 



H 



Similarly one can obtain the chemical potentials for D < 
0. The procedure is completely analogous to that used 
for magnetic impurities (degenerate Anderson model in 
the C/ — > oo limit "with Zeeman and crystalline field 
splittings)^ 



The phase diagram for D = 1/3 and D = —1/3 is 
shown in Figs. 4 and 5, respectively. For these pa- 
rameters the level crossings are ai H — I and H — 2. 
The level crossings stabilize the I-f-III+IV and I+II+IV 
(II+III+IV) mixed phases over the I+II+III+IV mixed 
phase. As discussed above, the phase I+III+IV for 
H > 1 involves different condensates than for H < 1. 
All phase boundaries are the consequence of one of the 
four rapidity bands getting empty and, hence, a tran- 
sition involves a one-dimensional van Hove singularity 
"with the corresponding consequences on the density of 
states and lo^w-T specific heat. For small magnetic fields 
the phase is a mixture of t-wo-particle and four-particle 
bound states. For larger fields H, the four-particle bound 
states are only favorable if the density of particles is high 
enough. Also for intermediate magnetic fields the phase 
diagram ioi D > is very different from that of 13 < 0. 



V. CONCLUSIONS 

We studied an ultracold gas of fermionic atoms -with 
an attractive contact potential by solving the correspond- 
ing Bethe ansatz equations. We obtained the phase di- 
agram for a 5* = 3/2 in a magnetic field (/i vs. H) 
"within the grand canonical ensemble. Four elementary 
states can occur: (i) polarized unbound atoms "with spin- 
component Sz = 3/2, (ii) bound pairs of atoms "with spin- 
components Sz — 3/2 and Sz — 1/2, (iii) bound states of 
three particles with spin-components Sz — 3/2, Sz — 1/2 
and Sz — —1/2, and (iv) bound states of four particles, 
one with each spin-component. Mixed phases of different 
classes of bound states dominate the phase diagram. For 
a given chemical potential the phases are homogeneous 
and display no long-range order. The transitions between 
phases are crossovers of the Prokovskii-Talapov type. 

There are several advantages of working in the grand- 
canonical ensemble vs. the canonical ensemble,~ where 
the number of particles is kept fixed. (1) By rescaling 
all quantities in the integral equations for the dressed 
energy potentials, e^"^', with the interaction strength |c|, 
one obtains universal equations for the phase diagram /i 
vs. H . Our phase diagram shown in Figs. 1 and 2 is 
then valid for all attractive |c|. This is not the case if the 
total number of particles is kept fixed. (2) Since the di- 
ameter of the tube gradually changes with position from 
the center of the trap to its boundaries, the effective local 
chemical potential varies along the tube. Within the lo- 
cal density approximation this change can be represented 
by a harmonic potential and as a consequence of the x- 
dependence of /i there is an inherent tendency of phase 
separation,^ i.e. the trap is inhomogeneous. At different 
positions of the trap then different phases may be realized 
and a sequence of transitions should be observed along 
the trap. (3) Josephson tunneling between tubes and in- 
teractions between particles in different tubes,— i^ may 
give rise to a dimensional crossover from one-dimension 
to a higher dimension. This gives rise to long-range or- 
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der of quantities that are generalizations of Cooper pairs 
for S = 1/2 to higher spin. The system remains strongly 
anisotropic and pure (there are almost no impurities) and 
is hence favorable for inhomogeneities like modulations of 
the order parameter of the FFLO type in the presence of 
an external magnetic field. All these properties can only 
be studied within the framework of the grand-canonical 
ensemble. 

We also investigated the theoretically interesting situ- 
ation of a quadrupolar splitting superimposed with the 
Zeeman field. In this case, the spin energy levels dis- 
play two crossovers as a function of the magnetic field. 
Hence, the character of the bound states before and af- 
ter the crossover changes. In general bound states of four 
particles are not favorable for low particle density or high 
magnetic fields. 

The phase diagram for larger spin values can be ob- 



tained the same way but the calculations and the results 
are considerably more involved and complicated. An ex- 
tension to hard-core bosons (particles with integer spin) 
is also possible^^i^ii^ since the Bethe ansatz equations 
for hard-core bosons are the same ones as for fermions. 
Note that wave functions for hard-core bosons have a 
node when two particles are at the same position, in anal- 
ogy to what is required by Pauli's principle for fermions. 
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